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THE BOUNDARY LAYERS IN FLUIDS 
WITH LITTLE FRICTION* 

By H. Blasius 


INTRODUCTION 


1. The vortices forming in flowing water "behind solid "bodies are 
not represented correctly "by the solution of the potential theory nor 
"by Helmholtz’s Jets. Potential theory is unable to satisfy the condi- 
tion that the water adheres at the wetted bodies, and its solutions of 
the fundamental hydrodynamic equations are at variance with the obser- 
vation that the flow separates from the body at a certain point and 
sends forth a highly turbulent boundary layer into the free flow. 
Helmholtz’s theory attempts to imitate the latter effect in such a way 
that it Joins two potential flows. Jet and still water, nonanalytical 
along a stream curve. The admissibility of this method is based cn 
the fact that, at zero pressure, which is to prevail at the cited 
stream curve, the connection of the fluid, and with it the effect of 
adjacent parts on each other, is canceled. In reality, however, the 
pressure at these boundaries Is definitely not zero, but can even be 
varied arbitrarily. Besides, Helmholtz’s theory with ItB potential 
flows does not satisfy the condition of adherence nor explain the 
origin of the vortices, for in all of these problems, the friction 
must be taken Into account on principle, according to the vortex 
theorem. 

When a cylinder is dipped into flowing water, for example, the 
flow corresponds, qualitatively, to the known potential, but as the 
water adheres to the cylinder, a boundary layer forms on the cylinder 
wall in which the velocity rises frcau zero at the wall to the value 
given by the potential flow. In this boundary layer, the friction 
plays an essential part because of the marked velocity difference; 
on it also depends the extent of the velocity-decreasing wall effect, 
which must be conveyed by shearing forces into the fluid, that Is, the 


*"Grenzschichten in Flussigkeiten mit kleiner Reibung." 

Zeltschrift fur Mathematik und Physik, Band 5 6 , Heft 1, 1908 , pp. 1 — 37* 
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thickening of the houndary layer. That the outer flow separates at a 
certain place, and that the water, set in violent rotation at the 
boundary, leads into the open, must be explainable from, the processes 
in the boundary layer. 

The exact treatment of this question was undertaken originally 
by Prandtl (Yerhandlungen des intern. Math. Kongress, 1904) . This 
explanation of the separation is repeated below. Since the integration 
of the hydrodynamic equations with friction is a too difficult 
problem, he assumed the internal friction as being small, but retained 
the condition of adherence at the boundary surface. In the present 
report, several problems, based upon the simplified hydrodynamic 
equations resulting from Prandtl’ s article, are worked out. They 
refer t'o the formation of boundary layers on solid bodies and the 
origin of separation of Jets from, these boundary layers suggested by 
Prandtl. The writer wishes to thank Prof. L. Prandtl for the sugges- 
tion of this article. 

2. The constant of the internal friction is assumed small aB in 
Prandtl’ s report. The boundary layers then became correspondingly 
thin; the fluid maintains its normal (potential) velocity up to 
near the boundary surface. Nevertheless, the decrease in velocity 
to value zero, and, as the calculation will show, the separation 
in this boundary layer must, naturally, continue, and so the potential 
flow is not completely regained, even at arbitrarily little friction; 
rather the separation and, the transformation of the flow effected 
through it behind the body must prevail even at arbitrarily small 
friction. 


Figure 1 


The procedure is limited to two-dimensional flow and coordinates 
parallel and at right angles to the boundary (arc length and normal 
distance). In spite of its curvature, the type of the basic equations 
in the narrow space of the boundary does not differ perceptibly from 
that for rectangular coordinates. With e as order of magnitude of 
the boundary— layer thickness 
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The friction gains influence when it is put at k ~ e ; this 
gives the relationship between boundary— layer thickness and sm a lln ess 
of friction constant. In the first equation, the term S^u/Sx 2 
cancels out; in the second equation, only Sp/Sy rsj £ or, when 
allowing for the coordinate curvature, ~ 1 remains . In both cases, 

^Allowance for the curvature of the coordinates produces, as is 
apparent when reforming the differential quotients, only in the second 
equation a not— to-be— neglected term pu^/r if r is the radius of 
curvature. This term is of the order of magnitude, unity. 


k 


me a tm 1256 


the effect of the pressure on y is to be disregarded since, in the 
narrow space of the boundary layer, the integration of $p/dy can, 
at the most, produce pressure differences of the order of magnitude 
€2 or s, or, in other words, pressure and pressure difference 
dp/chc are independent of y, hence, are "impressed" by the outer 
flow on the boundary layer. The velocity of the outer flow next to 
the boundary layer is denoted by u and. is to be regarded Bolely 
as function of x because the really existing dependence on y, when 
compared with the substantial variations in the boundary layer itself, 

can be ignored in the sense of the foregoing emissions; v is accord- 
ingly ~ e = 'fkj hence becomes zero with k. The remaining fundamental 
equations for the boundary layers are then: 



+ h. = 0 

Sx dy 


Boundary conditions are 

for y = 0: u=0 v = 0 

for y = »: u = n 

These equations establish, to a certain extent, a basis for a 
special mechanics of boundary layers, since the outer flow enters only 
in "impressed" manner. 

3. The qualitative explanation for the separation of flow 
according to Prandtl is as follows: the pressure difference, and 

with it the acceleration, is, apart from the friction term, constant 
throughout the boundary layer, but the velocity near the wall is 
lower. As a result, the velocity here drops sooner below the value 
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zero for pressure rise than outside., thus giving rise to return flow 
and jet formation, as indicated "by the velocity profiles in the figure 
below. 



Figure 2 


The region of separation itself is therefore characterized "by 


du 

dy 


= 0 


for y = 0 


This explanation does not work like the Helmholtz jet theory with 
an ad hoc assumption, hut only with the concepts forming the basis 
of the present hydrodynamic equations. The stream line, which bounds 
the separated part of the flow, departs at a certain angle frcm the 
area of separation since the stream function develops 

around the separation point [x] in the following manner: 


' t = Cjy3 + c 2 (x - [x] Jy 2 

As a less important effect, it is to be foreseen that, as a 
consequence of the stagnation of water effected by adhesion, the flow 
is pushed away frcm the body. Through this and the reformed flow aft 
of the body, the flow upstream frcm the body is, of course, affected 
also, so that the assumption of potential flow is Insufficient for 
quantitative accuracy of results and must be replaced by experimental 
recording of the pressure distribution. 


A. 
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I. BOUNDARY LAYER FOR THE STEADY MOTION ON A FIAT 
PLATE IMMERSED PARALLEL TO THE STREAM LINES 


The flow proceeds parallel to the x— axis . The plate starts in 
the origin of the coordinates and lies on the positive x-axis. 

In this very elementary case, there is no pressure, difference; 
hence, no separation is expected. However, the calculation is 
carried out to illustrate the mode of calculation to he used later. 
The fundamental equations read: 


P 





^+^ 1=0 

dx $y 

The equation of continuity is integrated hy introducing the stream 
function t : 


u 



v 


Sx 


Boundary conditions are: 


for 

0 

11 

£ 

11 

0 

v. 

0 

II 

> 

for 

y = °°: 

u = u. 

< — constant 


1 . According to the principle of mechanical similitude, the 
equations can he simplified when a similitude transformation converting 
differential equations and boundary conditions are known: multi- 

plying x, y, u, v, t hy the factors Xq, y D , u 0 , v 0 , and to results 

in 


P^o _ t _ u^o 

= y c ^ V ° = *> ; 


to = UoyoJ Uo = u 


as conditions that the problem and its solution are transformed, and 
that, through the trans format ion, p, k, u = 1 are created. The four 
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equations still leave a degree of freedcm in the choice of the 
factors x Q , y Q , u Q , v 0 , and t Q . The last three equations define 

the factors assumed by u, v, and iff through the transformation; the 
first states that the desired solution of the problem transforms in 
itself , provided only that 


pu £0 
k x 0 


= 1 


or in other words, with consideration of the factors which u, v, 
and \jr assume, the condition can depend only on 

pu y£ 
k x 

By this argument, the number of independent variables is reduced. 
Next 



are introduced; £ is then sole function of | and 


u = 1/2 

Insertion in the differential equation gives 

= £«l» 


Boundary conditions: 


for 

I = 0: 

g* = 0 

t = 0 

frcm 

£ 

11 

0 

for 

i = °°: 

CVJ 

II 

•• 


from 

£ 

11 

sh 


”44 
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2. The integration of these and subsequent equations is effected 

by expansion in series: expansion in powers for 5=0, asymptotic 

approximations for I = 00. The boundary conditions at both points 
being given, one and two integration constants, respectively, occur 
in the expansions. They are defined by the fact that both expansions 
must agree, at an arbitrary point. in the function value £, to the 
first and second differential quotient. The agreement of all differenr- 
tial quotients is then assured by the differential equation. 

3. Solution of the above equation by expansion in powers 

sr = -£** 

for 5=0 with the boundary conditions at thiB point 

V = 0 £ = 0 


is effected by 


t = ir ( 1 > n c n a ' n+1 t 3n+2 

5 & ( 3 n + 2 )! 1 

which is so chosen that the coefficients c n to be defined are 
whole positive numbers, which simplifies calculation. The 
n+1 

factor a brings out the nature of entry of the integration 
constant; c 0 , which otherwise would occur as such, can then be 
put as c 0 = 1. The recursion formula for c n reads 


'fife) CvC ^ v 

The first of the thus computed coefficients are: 

cq = 1 c-^ = 1 °2 = c 3 = 375 = 27,897 

C 5 = 3,817,137 2 Cg = 865,874,115 Orj = 298,013,289,795 

On account of the convergence, the denominator. (3n + 2)1 was used in 
the previous equations; £* and £" are easily formed. 

0 

The coefficients eg and c^ in the original thesis are 
incorrect. This error has no effect until the fourth decimal. 
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4. There is an additive integration constant for £ in the 
asymptotic approximation of g "because 

for | = co ; £ 1 = 2 

hence , 

£ = 2| + const. = 2r) 

so that T) appears as new coordinate shifted toward g . 

To compute a first correction £j_, put 

5 = 2ti + £ x 

which gives 

anti" - -£3.* 

with the squares of the corrections disregarded, hence "by integration: 



The general procedure for computing the other terms is such 
that further minor corrections £ n are added and its squares dis- 
regarded. The result is a set of linear differential equations for £ n , 

the left, homogeneous side always the same; at the right, the error 
appears as "impressed force" which the sum of the preceding approxima- 
tions, inserted in the differential equations, leaves. 

5 . The object is reached much quicker "by the following argument: 

The differential equation for £^ 

2^” = -£i m 

arises from, the original equation 


esT - -t" 
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when the roughest approximation £ = 2 tj is Inserted at the left for £ . 
Obviously, £ has the least effect at this point, and the differential 
equation is then integrated as if £ were hnown at this point. 

£ . / VV-JVn 

The three integration constants are contained in the arbitrary law- 
limits. Putting £ = 2 i\ at the right gives £-|_ at the left, but 

putting £ = 2 tj + £q at the right gives 

£ =j\r i y ,T 1 dT 1 e- T l 2 e-/ 1 ? l dT i 

or with consideration to the boundary conditions 
£ = 2q +7 f\ y^dqe-^l 

= 2 t] + fc _ 7 At, /V^r^dn 

i/00 £/oo C/00 


Hence, the second asymptotic approximation 

£2 - -7 2 /\ /’an . e"^ 2 /'’an /'’an /V^an 

Joo (loo Uoo (Joe Uoo 

By partial Integration 
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6. A general statement at out such, integrations reads as follows: 
According to the formula 

f e“ r l 2 T] n dT] = — iri n— ^e -1 ! 2 + — ~ P e~~ r ^r\ Dr -^d.r\ 

C/e 0 2 2 < ' c0 

to he gained by partial integration, each integral of this form can he 

__2 _2 

reduced to the functions e 1 and / e 1 dr) multiplied ty powers 

of t] • After several such integrals are obtained, the innermost is 2 
transformed, if necessary, in the indicated manner. The integral e — ^ 


AH 2 

or / e“ ^ dri multiplied ty powers, appears then telow the penultimate 

C/co 

integral sign. The former gives no new difficulty; the latter can te 
reduced ty partial integration tp the two functions e ^ and J e — 1 dr} 

/ 1T 'e _T l 2 dTi = t \ f T1 e -T l 2 dTi + ^O -1 ! 2 
• L (/« Uca 2 

J'Vn pe-^d-r] = iqS J T 'e- T l 2 dTi - £ pe~^ 2 d n + i 

pn^dT] p e~ 1 l 2 dTi = ^3 p e -T l 2 dri + ir^e-^ 2 + |e _Tl2 


ri e 


and so forth. 


If, as above, the integral can te quadratic in e~ ' H , four types 
must be distinguished: 




.-T / e -^ dr,, 


e ^ dr[ ' 


71 -2n 2 

dq 


multiplied ty powers of t} . The first and fourth types give nothing 
new. Partial integration provides for the second the formula 

fv^d, /V) 2 dn= j/Vi 2 d,} 2 - /Vi 2 !,, 

Uoo (ioo C/co 1 I WOO uoo 
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or 


V^dri pe~^\ = m'e-'W, 


and 


T]e 


-T) 


11 - 2 dn [V^d, =-ie ^ 2 Pe-^dt! + i 

J 00 2 Uoo 2 000 


2 e- T l 2 dn Pe-^dri = - ^ne^l 2 P e"^ + f J Pe^d^ > - Je" 2T l 


8 


and bo forth. 

Likewise for the third type 


,T,j 'VfiJ"*, - ,1 


Joo Woo 


!r-*4 


"IJ /"> -*2. 1 2 , 1 


700 uoo 


e _T1 dTjf T)dT] = 


/W 2 + |rie-^ 2 /V»% 


."¥1 


TJ [ 2 

e^ 2 dTi J- + ie -2 ^ 2 and so forth. 
0 


Since no new types for integrals are introduced ty these formulas, 

the indicated tables of formulas govern all integrals in which e”^ 
occurs no more than twice. Any number of successive integrations 
over such functions are possible; the powers of t) involved are 
unrestricted. The formulas for £ 2 3- n section 5 were obtained by 

this method. These integrations will be met again later r With the 
type of integration results thus known, the calculations can be made 
by utilizing a formula with indeterminant coefficients. 

7. With this differential equation, it is possible also to define 
the error that afflicts the present solution as a result of the effected 
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emissions . It is easily verified that 2r\j 2 t\ + tj, 2 ti + S x + £ 2 
remain 'below the true value of 5* An upper limit can also "be found 
"by employing the previously given (see Section 5) form.., somewhat 
modified 


£ = 2t) + 7 


/\ / V -^ 


( £ -2T])dTi 


of computing a finer frem a rougher approximation: a rather arbitrarily 

chosen upper limit, such as the first term of the semiconvergent 
expansion of for instance, is entered for £ — 2t), thus 


,-T 




and an asymptotically finer upper limit for £ is computed 

from this assumption. It is insured so long as the latter remains 
below the assumed one. The calculation gives (according to the general 
formula) 



4.1 £f. 

2 „V+1 



- 21^ < i fi-i 

0 n-3 


= ■d 


-/ 


(5 -2t)dTl 


< e 


e 
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Figure 3 


The upper limit for e^ is then found according to the above figure as 


e 13 1 + cfd 


g = 



1 


where t3q la the highest existing value of d, hence corresponds to 
the value of the coordinates for which 5 is computed. It 

results in 


£ < 2r\ + 7 



gy e -T 1 g \ 
8 n3 / 


< 2n 



ay 2 q " 211 

m n 5 
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similarly for £ * and : 


5* > 2 + $1- 


ay 2 e- 2 ^ 2 
32 r\h 


C" < Ci" 


+ 


cry 2 e - gl ') 2 

® n 3 


A more accurate execution of the integrals affords a more accurate 
result . 


8. The connection of the two developments and the determination 
of the integration constants (a, y s and r| — 5) is as follows: To 

separate the integration constant a. 



is introduced in the power development (3), which results in 


= =Y~(- l) n ^ j 

^ U ( 3 n + 2 ) t 


■3n+2 


dZ 1 t* - S / 1 °n v 3n+l 

= ^ (3* +1)1 


n=0 


(3n + 1)' 


= U" = Vj- l) n 7 ^ VT X 3n 


cLX‘ 


n=0 


(3^) ! 


The displacement of tj relative to g is expressed "by introducing 
the integration constant (3 


N^aT-Ti = X - 3 
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The formulas are completed "by inserting 

£= 2tj + £-]_ + £ 2 

5 ' - 2 + + Z 2 ' 

£" - ^ + £ 2 " 

from the asymptotic approximation ( 4 ) and (5). A graph is made for Z 
and its differential quotients from which the following values are 
quoted: 


X = 

0 

0.8 

1.0 

n 

m 

EH 

2.0 

2.05 

mm 

z = 

0 

0.317 

0.492 

0.701 

0.938 

1-63 

1.79 

1.8561 

1.94 

dZ _ 
dX 

0 

.784 

.961 

1.121 

1.257 

1.50 

1-53 

1.5479 

1.56 

= 

dX 2 

1 

_ 

— 

— 

.639 

• 3 ^ 

.28 

.2582 

•23 


The terms of the power series are computed up to 


Zr jX.‘ 


■23 


23! 


further terms 


are extrapolated^, in part^ from the difference series of the logarithms 
of the coefficients. The location of the asymptotes is already quite 
apparent in figure 4 . 
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Since owing to £ = 2t] asymptotic. 


a2/3 


•(X — P), rough approximation 


values can already he read for a and P: a = 1-30, p = O.96, 

with X = 2.O5 as connecting coordinate for rj = 1.00. The corresponding 


values of and give for 7: 7 = O.92. The calculation is 

dX 2 

more rigorous when a ,7, and the connecting coordinate related 
to tj — 1 are varied hy minor corrections and these then computed 
from linear equations. T p judge the accuracy, it is stated that our 
calculation for X = 2.05 gave 


z = 1.8561, §§■ = 1.5^79, = 0.2582, = -0.479 

“ dX 2 dX^ 

where the fourth decimal is no longer certain. Asymptotic approximation 
for 7=1 giveB (using Markoff’s J e — fc2 dt) 
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£ =2 + 0.0l*;*-7 +f 0 :™^7 2 

5* - 2 - 0-139W • r - {S:X 7 3 6 } r 2 
5" = 0.36788- 7 + |°;^|} r a 

the top numerals in the | jstemming from £qj the bottom numerals from 

the upper limit defined in (7)* (The latter is, as stated "before, 
rather rough.) The "temporary assumption" about the upper limi t of £ 
gives: £ < 2 + O.O92 y. The upper limits are therefore guaranteed 

(reference 7) . Hence, the result 

a = I.3266, X = 2. Okgk, y = O.9227, (0 = O.9508) 

It can be safely assumed that a ranges between I.326 and 1.327. 

9. From it, it can be computed, for example, what drag a plate of 
width b and length 2 is subjected to when dipped parallel to the 
flow lines into a flow moving at velocity TU The drag per unit of 
surface is 


l = ,m -i 


y 


dy 


VE 



JhpP" -jL 


Integration over the plate gives 


b 


IS 


= Sib 
2 



hence, when the water flows at both sides of the plate 
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II. CALCULATION OF REGION OF SEPARATION BEHIND 
A BOB! DIPPED INTO A UNIFORM FLOW 


1 . The following problem Is treated: In an otherwise parallel 

flow, a cylindrical body Is Immersed symmetrically to the direction 
of flow. The boundary— layer coordinates are computed from the point 
of division of the flow. The quantity u is expanded as function 
of x in a power series. For the integration of the fundamental 
equations 


..du , _du _ -rdu 

U-«r— + Vrr — = Usr— 

OX OJ ox 


, 'k S^u 

P Sy 2 




* - 21 izx 2Z+1 

l=o 

the formula 

l=o 

is used, with due regards to the symmetrical conditions for the 
stream function u and v are obtained then by differentiation. 
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Consistent with the general "boundary conditions, the functions \(j) 
must then satisfy the boundary conditions 


hence 


Xj' = 0 X 2 = 0 for y = 0 
X * = q. 2 for y - 00 

■ v * r i 


.Vi is the constant of integration. From Insertion in the first 
fundamental equation, the differential equations for X are obtained 
as: 


+ IXYW - (a. + + #V 


1=0 


1=0 


which for 2=0 is quadratic, for 2 > 0 linear in the X 2 function 

to be' defined. This equation can, like the proceeding problem, be 
solved by expanding y = 0 in powers, for y = « approximating 
asymptotically and Joining both. Subsequently, it is shown that the 
asymptotic approximation can be omitted, since the power series already 
identifies the asymptote and therefore the integration constant with 
sufficient accuracy. The calculation is restricted to x q and X]_> 
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that Is, the first and third powers of x. Because, since the corre- 
sponding coefficients q Q and ^ In u already indicate a first 

increasing, then decreasing velocity — the case, in which presumably 
separation occurs, is characterized ty q Q > 0 , < 0 — the type 

of pressure distribution required in the introduction ( 3 ) is already 
supplied "by the first two powers; hence, it is to he expected 
that Xq and even though not quantitatively exact, already 

represent the effect of the separation. Xn one of the problems 
treated in similar manne r later on the next approximation was also 
computed; and it substantiated the admissibility of the limitation 
to the first two powers of x. 

2 . The equations for Xq ^ are 


V 2 “ *0X0" = lo 2 + jyXo 



*o*L 


+ 3(X 1 *X 0 T - XjX") = 4q o q x + 


*•* 

P X 1 


The manner of entry of q Q , k, P can be established by mechanical 

similarity. Here also, the first two terms indicate universal 
significance in same respects. Hence, writing 


u = 


V 


*l x “ 


t = 


X x 

o 



and introducing the following quantities 



for x, y, x D , Xq gives 
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* =l/fe 5 ° 51t i 53 > 


^ 1 ES t * 


u ^ ~ 


2«a n 


( C c * I ± C^* |3) etc. 


£ o and ^ satisfy, as functions of T|, the differential equations 


Co 2 -C0C0" - * - Co" 


«oV - 3 £o” 5 i - S 0 Sj" - 16 + Si 


I* I 


Boamdary- conditions 


for 

n - 

0: 

Co = 0 

So' 

= 0 

O 

t! 

1— 1 

Xj) 

Si' 

- 0 

for 

n = 

00* 


So’ 

= 2 


Si’ 

= 2 

3. For 

So 

the 

power series 







c = 


a^ +1 b, 


ta=2 




is entered. 

Insertion in the differential equation gives: 

l>2 arbitrarily = 1, since a already is integration constant. 

0,^(3 3 = —4; since, in the formula of the integration constant a, 
no allowance was made for the homogeneity of the equation fern* £ 0 , 
a appears again in this equation. 

b^ = 0; the curvature of the velocity profile does not change, 
at first, since the friction in its effect is two terms ahead of the 
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inertia; starting frcm m =.5, it is 


m 





"bp."bm — — p. 


The coefficients of these recurrence formulas can, like all numbers 
combined this wa j frcm binomial coefficients, be computed from, a 
diagram similar to Pascal’s triangle, whose start is the following: 



and in which each term is the sum of those above it. Only the 
framed— in portion, consistent with the foregoing limits of sums, 
is counted. 

The first 13 coefficients are 


b 2 - 1 

■*\% 

1 

11 

on 

& 

\ = 0 

b 5 = ! 

*6 “ 2b 3 

*7 ’ n j* 

bg = -1 

”9 = -*3 

b i0 = -l6b 3 2 

b~i i = 27 — 16b g ^ 

b 2 2 = lBlb 3 

b l3 = 840 b 2 

Besides a, two more integration constants due to the 
asymptotic approximation are involved, which, as in the preceding 
problem, should join the ccnqrated power series. For the present 
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purposes (calculation of point of separation), it is, however, 
sufficient to know a, and, as stated before, it will be seen that a 
can already be computed with sufficient accuracy by means of the 
power series. 

1 cLZ 

Put Z Q as — Cq, H = arj and plot — as function of H from, 
a cLU 

dZ 1 

the power series. — —■ itself is still dependent on bo = — — and 

dH 0 an 

shall, for the correct value of a, approach the asymptote . 

dH a 2 

For other values of a, it approaches no asymptote at all, as a result 
of which as fig. 6 shows, the method for defining a is very 
sensitive. The value a = I.515 is obtained; the last cipher is no 
longer certain. 


dZ 

o 



5. The calculation of £ "by the above linear equation and the 
boundary conditions is effected in similar manner: power form ula 
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C2 = 1, since 8 already .is integration constant; 8c 3 = — l6; 
cj. = 0; and for m > 5 : 



Here also the coefficients in these formulas can he computed 
from a diagram whose first line (m = 3) consists of the numbers — 1, 
+4-, — 3, while the others follow by addition: 



The first coefficients are 

c 2 = 1; 8 c 3 = — 16; cj^ = 0; c^ = 4a 3 ; eg = 6a 3 c 3 - 8; 

C’j rn -32c 3 J eg = 17a 6 ; °9 = 30a 6 c 3 - 224a3; 

c 10 = -576a 3 c 3 - 25 6; c^L = 2048 c 3 + 294a 9 ; 

Cj £ = 783a 9 c 3 - 5092a 6 ; c l3 = -17392a 6 c 3 + 59648a 3 ; 

c x ^ = 221952a c 3 - 315a - 136192; 

a = _11025a 12 c o - 1024000c, - 54864a 9 ; 

15 3 3 

c l6 = 174l68a 9 c 3 - 221296a 6 . 

6. The asymptotic approach is again disregarded, the integration 
c ons tant 8 being defined by the condition that £3/ must have the 

asymptote ^ =2. Figure 7 shows the terms of the power series 
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for £ ]_* , those free from and those multiplied by eg, as - 

curves A and B, that is 

= 5 • A _ 16 - B 

after which is plotted for different values of 5. 


£ ' 



This curve indicates that the convergence of the series is rather poor 
in spite of the great number of computed coefficients c, even 
at T] = 1.6. In any event, the terms indicate, when identical powers 
of a are combined, a satisfactory variation so that the series are 
still practicable. The correct value of 5 ranges between 8.20 
and 8.3O. The curve rises, at first, very quickly and approaches 
its asymptote from above. This marked influence on u near q = 0 
compared to t\ = <» permits u in the case of separation to change 
signs at the boundary before it does on the outside. 

7. Proceeding to the calculation of the point of separation, it 
will be remembered frcm (l) that, quantitatively, the results are not 
exact, since only the first and third powers of x were taken into 
consideration. The point of separation [ij is defined by 
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or t>y (3) and- (5) 

a3 ± 5 [l] 2 = 0 

By (4) and (6), respectively, a = I.515, 5 = 8.25- Thus, in the case 
of the lover prefix, the only one of interest, the coordinate of the 
point of separation is 

[|] = O.65 

hence 




with 


u * q. 0 x - q x x 3 


The maximum, of the velocity (minimum pressure) lies therefore at 


x 


0-577 



while zero velocity in the outside flow would not "be reached 


till x * 1 • \f— • Accordingly, the point of separation is 12 percent 
a l 

of the total "boundary— layer length "behind the pressure maximum. The 
obtained figures are independent of friction constant, density, and a 
proportional increase of all velocities. 


According i to Prandtl* s diagram (section 3 of Introduction) the stream 
line = 0 diverges from the "boundary at a certain angle, which is 
computed as follows: In the vicinity of the point of separation, the 

development of the expression for if given in (2) reads 
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Ho 1 2 1 


♦- Usir F ((5 o m H - 5i” [i] 3)t i 3 + 3(Co” - 36,- S 2 )(I - [I] h 2 ) 


^ = 0 gives for the divergent stream, line 

2 3 

— rr. J 


= 3' 


35 


a 


i-.jj] 16 [i] 3 - ^ OQ 


= 11.5 


or in the not— reduced coordinates 


x 



H.5 



These formulas are characterized by considerable uncertainty because 
only two terms of the development of were computed and the higher 
differential quotients, which represent more subtle processes, are 
always less accurately computed than the former. 


III. FORMATION OF THE BOUNDARY LAYER AND OF THE ZONE 
OF SEPARATION AT SUDDEN START OF MOTION FROM REST 


1. The two preceding problems treated stationary flows. The 

problem of the growth of the boundary layer is now treated. Assume 
that a cylinder of arbitrary cross section is suddenly set in motion 

in a fluid at rest and from t = 0 is permanently maintained at 
constant velocity. At first, the state of potential flow is reached 
under the single action of the pressure distribution. The thickness 
of the boundary layer is zero to begin with, so far as the sudden 
velocity distribution can be obtained at all. The boundary layer 
develops in the first place under the effect of friction, then 
through the convective terms. The result is that, after a certain 
time, the separation startB at the rear of the body and, from there, 
progresses gradually. Since the fundamental equations refer only to 
thin boundary layers, they, naturally, represent only the start of the 
separation process, Just as the previous problems dealt with the boundary 
layer only as far as the zone of separation. 
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2. The equations involved here are 


h + u|Si 

dt ox 


+ v- 


chi 



+ K 


^U 

Sy 2 


= a* = _ 

dy ox 


K 


lr* 

substitutes for 


The potential flow which is set up first gives the houndary 
value u as function of x. Since the process for t = 0 is singular, 
the type of development is, for the time being, still unknown; it 
must be established by successive appro xim ation. The principal influence 
on the changes has (at small t) the friction, hence, for the first 
approximation u 0 


du 0 

at 



dy 2 


The integral of this equation 


u 


o 


_ 2u 

-JS 



dq 


*1 = 



satisfies the conditions of supplying a vanishing boundary layer 
for t = 0 and of joining the outside flow u G = u for y = 00 . The 
subsequent approximation is obtained by Inserting u Q in the 
convective terms, while time and friction terms obtain u = u c + Uj_. 

The resultant equation for u^_ reads 


c^i aV 



(function of q) 
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According to mechanical similarity, this equation is satisfied by the 
formula 


U 1 - f(n> 

which is also not contradictory to the boundary condition u^ = 0 
for y = 0 and y = ». 

After further considerations, vhich in particular refer to the 
insertion of x, the quantity u is represented in an expansion in 
powers of t, the coefficients of vhich are functions of tj, that is, 
still contain t. These functions are also still dependent on x, 
but this time x enters the differential equations only as parameter. 

3« The formula for ifr is accordingly 


t « 2^Kt Z^tVx^xi}) 

V=o 


= 


2^~Kt 


8*V 


U = H tV S 

VaO °11 


and hence the differential equations for X 






for (i = l, the right-hand side contains 

As before, the calculation is limited to the first two terms, 
that is 


X c - «e o (n) 


*1 - 



MCA TM 1256 


31 


hence 


u = u5 Q « + 

The equations for £ and . read then 

5 0 ” + 2t]? 0 " - 0 

? 1 "' + 2n£j" - ‘•? 1 ' - MS 0 ' 2 - £ 0 £ 0 " - 1) 

Boundary conditions 

for ti = 0: 5=0 5*=0 

1 3 o 3 o 

k - 0 k’ - 0 

for r\ = eoj S 0 * =1 = 0 


4. The solutions of the above differential equations, which are 
to he used in the subsequent problem, are obtained by quadrature when 
the homogeneous equations are integrated. The latter integrals were 
obtained by the following consideration: The homogeneous parts of 

the equations stem from the time and friction term which together 
form the heat conduction equation 

5u _ 

dt ” By 2 

Of this equation integrals of the form 


Un = t^Cn) 
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exist, according to similarity considerations, whereby f n satisfies 
the differential equation 


f n + 2 t iV " ^n = 0 
which the above form possesses. 

Thus, for example (see above) 



For t] = 0 

u Q * 1 when t > 0 
u c = 0 when t < 0 

For Tj » 0, hence, for y = 0, is proportional to t n , hence must 

be representable by superposition of solutions u 0 in the following 
form 


"n “ n 


>00 

u, 


= n u 


\2 ^c(t - t 0 ) / 
^k(t - t Q )) 


* 0^0 


to^ito 


10 \2\flct 


since for y = 0 it is 




dt 0 a t 11 


For the evaluation of this integral, put t — t 0 = t 
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insert herein. 


2\pct 


2 \j kT 


and finally obtain 


4- _ 1 jfi. _ - 

t " W ~ 2 S T " IS 72 * 


d.T = — — ^<34 ; 

2k £3*' 


Calculation of this integral by the binomial theorem and the 
previously cited method of partial integration finally gives 


n 2^ 


fn - 2 —Tz\ t 

u=o v ^ 


n n 


V ^ M 

1 J... 3 . 1 1 (/» 


0 


+ JZ<T ^ >t 1 2v_1 ®” ti2 

v3 ]t^ (2n - 1) . . . (2n - 2v + 1) 

The other integral is algebraic and equal to the above factor 

P n _o 

of / e~ ^ dT] 


f n 2 ^ 

n f£ 5 ( 2 n - U...3.1 


5. Quantity is determined as follows: With the boundary 


conditions taken into consideration 
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■whence "by Integration 


t = - -pr + T] + 
b o fH 



lT1 e-^ 2 dTi + ±e^ 2 ' 
2 


while utilizing 


p " _ 2 0 ~^1 £ 

“f? 


The second differential equation (of the second order for £ ^ ) 
assumes then the form 

I ~( \ r 

u r.^>* 

it 


16 


V> 2 a„ - Ane - 112 ♦ ±£ 
V 


{/.'■'•'I 


- T|6-^ 2 JV^di, - ie-^l 2 + £rv 

The integral of the homogeneous equation for £ 1 is "by (4) 


f * a(2r) + 1) + P 


r^e - ^ + (2^ l) e - ^ dr) 

II CO 


The Integral of the nonhamogeneous equation would then he 
obtainable by quadratures. But it is also true that, by twice 
differentiating, the differential equation finally becomes 


£ Mn * + 2t]^ ,,n = function of Tj 



mCA tm 1256 


35 


whieh IB easier to integrate as an equation at esseatiaW «** 
order . Hence 


{; nn = e“' r T 


/^e*! 2 [function of T) ] OTj 


Since tta pressed t~ce ot tte di«erential equation contains e- 

...nj-u cancels out, and t 

SrtSnT c^Ste^tc^te inactions teMni tte 

e"^ 2 twice, and in addition, powers 
integrals contain, several times, which can he accanr- 

of -n, and must he integrate - vy. (j 6) . The result of 

■nlishad hy the methods discussed previously * 

S StteJ roluminous calculation reads 


h' 


* J V^ 4 l + |' 2t|2 “ 11 


09 






n -ft 

2„„- £h l 2 


S 1 




_ i( 2 n 2 + 3 )s r ‘ S t ^> e 11 
\lic 


+ 1+ar} + P 


2 

2e~^ + 

[1 ^ „2 
© —T 1 dri 

L 

00 
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due to the fact that ^ is simpler than u~i, although "both have, 

at ax 

according to the order of the differential quotient, "heat conduction 
character." 


The determination of a and 0 frcm the limiting 
conditions f^_* = 0 for tj = 0 and t ] = » gives 


a = 0 



+ 


k 

3 ^ 72 * 


6. For computing the zone of separation, there is 


2 \|rct \ 0 ° x 1 / 


far ti = 0 


Then 


t " - 2_ 
0 " 


.< _ 2 


8 


?1 ” * ^ " 5572 


The condition for the time of separation [t] is 


1 + 




hence, ^ must he negative. The separation occurs first where ~ 
dx ox 

has the greatest magnitude. The result applies to cylinders of any 
cross section; u is the corresponding potential flow. 


I V. DEVELOPMENT OF ZONE OF SEPARATION FROM REST 
AT UNIFORMLY ACCELERATED MOTION 


1. Against the physical principles of the foregoing problem, the 
objection may be raised that the sudden shock might be accompanied 
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"by an Interruption of the fluid. Hence, let the solution of the 
problem assume that, starting frcm the time t = 0, the immersed 
body is subjected to constant acceleration. 

In that case 


u = tw(x) 


1 dp 3u _Su , 2 

+ u_ s W + t c W— - 

p ox at ox ox 


2 . From considerations similar to those made before, the solution 
of the differential equation 


^ = - i & + 

St Sx Sy p Sx ^^2 


is based on the formula 


ljr = 2 . ZZ' fc2v+:L>c 2v+i( :x: n) 

v=o 


CO 

u = ) t 


v=o 


2V+1 ^-2v+l 
Sti 


n--* 


\J0t 


Insertion in the basic equation gives 


S 3 X 


2A+1 


S 2 X 




+ 2tj— 


2X+1 




- 4 ( 2 X + 


)_ ^r 
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d^X-Sia-l _ ^ X 2^+l 

Stj dxSy Sx ^2 

for \ = 0, the right-hand side contains —I kw, for \ - 1, -4v£s. 

dx 

The calculation of the state is again l imi ted to the first two 
terms, while it should "be noted that through those two terms, the two 

terms of the pressure w + t^-w^S are also taken into consideration. 

ox 

The impressed force of the next equations contains only earlier 
development coefficients. For the final equation, however, which 
supplies the zone of separation, the coefficient of the next term is 
computed also. Fpr ^ and x^ the relationship of x can he 

introduced in the following manner: 



\ = w Sl (T l)', X 3 “ 


The differential equations for £ are then: 


&■, 

± + 2ti ± - Ui = -4 

ar]3 


Stj* 


+ 2q 1 - 12^—3 = -4 + h 






c ^\2 

Sri 


- C 


ti 1 
1 *1®. 


Boundary conditions: 
for q = 0: 


= 0 , 

^3 

as, 

for tj = «« — i = 1, 

dr\ 


£.0 

on 

Uo 

8^3 

Sti 


> 


= 0 


u = 0 


from 


v = 0 

from u = tw 
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3 . According to the general solutions of tlie present type of 

as, 

differential equations discussed in III (4), rr— = can be written 
forthwith, since the nonhomogene ous term. — k is disposed of by 


^n 


= 1 


is obtained hy integration by the repeatedly cited method ( 16 ) 


ftl.jL 

St] 2 fit 


2 pn 2 

e n +2t\ e“n dtj 


as 


l - 1+ l 


2 nn o 

t^e - 71 + (1 + 2t\ 2 ) e~^ dq 

U CO 


■ n + 


3^L 


- 1 + (1 + Ti 2 )e-^ 


( 3 n + 2 tj 3 ) J* Vn £ "dti 


These functions are quantitatively plotted in figure 8 and given in a 
table (see Section 6 following) . 



l+o 


NACA TM 1256 



Figure 8 


The impressed force on the right-hand side of the second equation 
is then 


16 

fit ^ 


P T 'e -T l 2 dri + 


16 

3* 


0^2 2 
4t] e 1 dt] + 2e ''l 

_ Uoo 


+ 


16 

3« 


(-2 + t] 2 )e -2T l 2 + ( — 4 t} + 1+T]^)e T l 2 


*1 »,2 

e ^ dT) 


+ (3 + W*) 



4. The integration of the second equation, in closed form, again 
succeeds by the same methods as in XII (5)’ For "k* 1 ® part of the 

2 

impressed force quadratic in e ^ a formula with indeterminate 
coefficients is particularly advisable. 
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if-1 


(a + t>T] 2 + cri^)e“^ T l 


_ c 2/^2 

+ (dr] + 6T^ + fT] :? )e ^ . e T 1 dtj 


+ (g + ht] 2 + Itj^ - + kt^) 


^ 2 , 
e ’l dtj. 


This formula falls when the Impressed force contains terms vhich exceed 
6 — 2 t|2 5 — ti 2 — 1 n^. k J P 1 ^ — n2. 

a 1 *n ''o 1 


T| 0 


t \ e. 


e 11 df], T]” 


e 1 dT) 


(compare III (5))* The 


coefficients are dete rmin ed frcm linear equations. The other 

S 2 ?-. 

portions of — — *■ are easier to compute ; £ 0 and 1 f ollov "by 

St 1 3 St] 2 

integration and differentiation. So, vhen the integration constants 
are correctly ccmputed, the final result is 




3 _ 


^ 0 — n 2 32 P 

3 ^r 1C5 « 


71 -vn2 

e“^ dTi 


3« 


Sti 2 3 jiT ^ ' J °° 

(-T) + 2T]3)e _2T l 

+ (1 + 2r[ 2 + 8r] il )e -T l 2 P e" J 'l 2 dii 

Uoo 

+ (6t) + 8tj3 + 8t^) 1 J e — ^ 


2 I 
dT 7 J 


+ ±p ^=\ , 

5 \fifn U5'Pt3 y L 


(16 + 3<Sn 2 + 8r 1 li )e-^ • 

Tl 


+ ( 60 t] + 80^3 + l6ri^) j e ^ 2 d.Ti 
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+ — (8 + T] 2 + 2T] 1; )e~ 2T1 

9n 

+ (2kr\ + 8t} 3 + 8ri^)e — ^ / e-^dTi 

V 00 


+ (-9 + 18 ti 2 + 12 t]^ + 81^) ) P © ^drj 


+ T^f~= T?) (33T1 + 28n3 + k ^ )e ~^ 

15 \ 6 ^ 45 ^/ L 

+ (15 + 9°t ) 2 + 60 V* + 8 ti^) J ^ e - ^ cLt] 


hence, "by integration 

t = - -2-Te-^ 2 + 2t, f \^ 2 d n 
3 3 'Iff - U Co _ 




Tje - ^ 2 + (1 + 2 t) 2 ) e -1 ! dr) 


(49ri + 11 t]3 + 10rp)e~^ r P + 768 P e” ST l 2 dTi 

Uqo 
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+ (-537 + 198ti 2 + 64^ + 40 n 6 )e^ 2 



+ C — 33-5n + 210 ti3 + &hr^ + 40tj^) 



105 45 P 


(24 + Qjrfi + 


40t]^ + 4Ti^)e — ^ 


+ ( 105 t) + 210ri^ + 84 tj^ + 8rj^) e -1 ! 2 !!) 


/ 128 128 _ 9 \ 

\!575^t3 105^2* 14\[7 j 


These three functions, plotted in figure 9 j rigorously satisfy 
the differential equations and the boundary conditions for the 
coefficient X. 


5 . The condition for the zone of separation has the form 




+ t% 

ox 



whence, by the foregoing formulas 



Si)2 yfj[ St] 2 15/rt 


256 

225^3 
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The equation for the separation time [tj reads 


1 + /ai _ [ t ] 2 |i£ = 

l 60 225 k J Bx 


The next term in the separation equation = 0 would 


read: — i— t^. 5., and in order to he able to allow for it, too, the 

2^t Bq 2 

coefficient in this separation equation, rather than the total 
variation of is computed. 



Figure 9 


The development term X^ satisfied the equation 

^3x 5 ^ B 2 )^ '^B^ ( ^ 2 x ± Bx x S 2 x 3 ^ ^ 2 x 3 Sx, 

^3 + 2 2 Bq Sq BxBq ox ^2 + Bq BxBq Bx ^2 

The entry of x in X^ and X^ is known, and calculation of the 
right-hand side confirms that X,_ assumes the form 
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X, 


■ -(!)%*•= S' 


5a 


Since tv cancels out, the condition of separation reads 




%T) C 


4- t 2 ^ 

Sx 


-lri=0 


a 2 s. 


fr -M> 


+ 1 4 )^ 2 

\3x 




r 


Jt ]=0 


+ t\£z 

cix 2 


5 2 6 


5a 




and 


6. This leaves the calculation of the coefficients 

~a 2 s ( 


Jt ]=0 

325 




152: 


&n s 


T1=0 


-It |=0 


For £5* the differential equation reads 


a% * 35 , 36, „ a 2 ^ , a?6, 

1 + 2t) S _ 20—5. „ &— -i —2. _ us . — 2 - 45, — i = 

*\ Q •x P r^n r^<n P 3 ^^2 


aq3 St ]2 a n 3i) 3 tj 


V 


£(t]) 


and the boundary conditions 

as 


S_ = 0, — - = 0 for ti = Oi 

5 an 


aSp 

an 


for tj => 00 


The impressed force f(rj) is given by the previously written 

a^sr 


functions. The desired coefficients 
Green* s method as follows: 


arf 


are computed by 


J tj=0 
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as,, 

+ 2 tv ^--2 _ 20—2 

V ^ 


dTi 



M 


+ 2tv3 


as 




J o 



Then, If -a is made to satisfy the adjunct differential equation 


6 2 6 

a n 2 


2t^ - 226 

' °n 


0 


and the boundary conditions 


6(0) as —1 6(°°) » 0 

the result is 

a 2 ?^ 



f ( iq) is given previously; the influence coefficient 6 (Green* s 
function) is obtained by integration of his differential equation 


6(n) = 


945 Jit 


(2895n + 5280^ 3 + 2352 ^ + 352n 7 + l6r\ 9 ) 


+ (945 + 9450tj 2 + 12600V*- + 5040ti 6 + T20ti 8 + 32 n 10 ) e n 2 PV^dr) 

QD 


The curve of 6 is shown in figure 10, along with the 
product 6 • f. The area of this last curve gives the desired 
coefficient 
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The computed values are the following: 


r\ 

0 

0.25 

0.50 

1.00 


00 

^1 

0 

.061 

.211 

.638 


n - °*376 

&Si 







an 

0 

.450 

.720 

.943 


1 

3rj 2 

2.26 

1.396 

•799 

.201 

O.O35 

0 

5 3 

0 


.060 

.115 

— 

.138 








i 

an 

0 

•137 

.150 

.020 

— 

0 

&i 3 

8t} 2 

.964 

.231 

-.092 

-.156 

-.05 

0 

f(n) 

0 

.315 

•750 

•457 

— 

0 

$(t]) 

-1 

-327 

-.112 

-.018 


0 

• f 

0 

-.103 

-.084 

-.008 


0 

e(n) 

0 

.124 

.240 

—0l8 


0 

•a • g 

0 

-.041 

-.027 

.0003 


0 


The area of the two curves is approximately 


'jfj 

%\ 

1 — 

« -0.058 


_ 

q=0 

Sq 2 

— 


T]*0 


-0.023 


7<> The equation of separation therefore reads 


i + [t]2k(jL - -256 ' ) - [t ]%) s 

{k 8x^15^ 225 ifl ->} \chc/ 


0.058 - [t] 


4 ..c£w 

chc 2 


0.023 = 0 
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or 


1 + 0.427 


{t] 2 |2 - 0.026 


HtiT - 


0.01 


Sx 2 


Since the newly added correction term is even negative* the existence 
of the zero position appears to he certain. 

The position and time of separation is apcording to the earlier 
approximation (without the term computed last) 



For the case of a c ylind er symmetrical to the direction 

flow* w^% » 0 at the rear point where the separation starts* the 
Sx 2 

newly computed correction gives 



* -2.08 


equivalent to an error of about 10 percent. From this the quality of 
the approximation made in the other problems* where only the first 
powers were taken* can probably be also appraised. 

V. APPLICATION OF THE EESUITS OF THE SEPARATION PROBLEM 
TO THE CIRCULAR CYLINDER 


1, On the circular cylinder 


u => 2Y sin § 

K 


£ is called the reduced coordinate X: V is the velocity at which 
R 

the parallel flow flows toward the right* and the cylinder moves 
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toward, the left., respectively. In the steady case, the separation 
starts according to Part II, Section 7 at Xg 0 p, = O.65 the 

maximum velocity lies at 




where 


u = q x — a x^ 
u o *1 



Taking the ordinary development in powers of sine 



'll 


27 

6r3 


sep. 


I.59 . B, X 


sep. 


- H 1 


"Siaj 


1.41 . E; X = 81° 
max 
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But approximating the slue in the Interval 0 — jt "by the method of 
least squares, gives 

<1, « ^ • 0.093 
x r3 


— . O.856 , 
r 


x 

sep. 


x 

max 


1-97-K, 


1.75 • R* 


X 

sep. 

X 

max 


H30 

101 ° 


In any case, the point of separation lies, by the present 
calculation, at from 11 percent to 12 percent of the total "boundary- 
layer length "behind the ma-HTmim of the velocity. TMb statement makes, 
of course, no claim, to accuracy, since only the first two powers of x 
are taken into consideration. Besides, test records of the pressure 
difference indicate that the state near the separation is difficult 
to attain "by development from starting point of the "boundary layer, 
Because it is too strongly affected "by the pressure distribution of 
the turbulent bodies behind the cylinder. The sole purpose of the 
present calculations is to indicate that separation is actually 
obtained by the hydrodynamic equations. Further development of the 
calculating methods, especially for the more important problems of 
solids of revolution, promises, therefore, success. 

2. If the cylinder with constant velocity is suddenly set in 
motion 


u » 2Y sin X 


5u 

dx 


27 

R 


cos 


X 


The time of separation 



is, according to III (6), given by 



[t] = -O.35 


R 


V cos X 
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The separation starts for X = it, cos X = — 1 at time 


0 . 35 ^ 


Up to that, the cylinder has travelled a distance 


S a Yt c »0.35 . B 


All this is independent of velocity, density, and friction coefficient 
(little friction assumed). 

3. At constant acceleration 


u b tw(x) 


2Yt sin — 
E 


where Y is then the acceleration of the cylinder in the flow. The 
separation time is (TV, 7) for the start of separation 

[t] 2 |^ = -2.3U or b - 2.08 


respectively, or 


ftp b - 1.17 S or 

Y cos X 


-l.Olt 2 

V cos X 


respectively. The distance covered hy the cylinder is 


s = lyt 2 

2 


at start of separation (X a it) 

S = O.59 . E or bO«52.E 


respectively. 
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4. Ttie resistance which, the cylinder experiences at constant 
acceleration is c omp uted next. The stress components are 

Y 7 - -» + ^ 


*7 



arifi ^31 


Owing to the smallness of the friction., ^ 

respect to ~ 3 leaving as force in direction of the 
qy 


dx 


cancel with 
outside flow 


K = 2 



cos X + 




RdX 


B is the width of the layer (height of immersed part of cylinder) . 
The pressure portion is computed as follows; 


■^pr 


assure 


TJt 

2BR p cos XdX 

do 



sin XdX 
ox 


Then 


dx 




u * tw 


p w + 



w = 27 sin X 
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The second term cancels out in. the integration; the first gives 


■^pressure “ 2 jtpBR^Y 


hence., an increase in inertia "by twice the amount of displaced fluido 
The friction portion is 


■^friction 



where k ■ k/p. Again, the second term disappears because 

and d are merely constants, leaving 

Sri 2 



^friction - 4 \f*pkt . BRV 


5 « To give a picture of the flow ccnditions corresponding to 
these formulas, the flow curves for a specific state of motion of the 
unif ormly accelerated cylinder are represented in a diagram. The 
parameters R, V, k are arbitrary; hence, necessitate the introduction 
of reduced quantities for x, y, t, and u, so that R, V, K 
disappear. It is accomplished by 
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■by which the formulas (compare IV (2) and. V(3)) 



■became the following reduced equations 


I- 



sin X 


( ^ + 2T 2 cos X 



TT « 21 sin X 



+ 2T 2 cos X 


*a\ 

) 



The curve $ is then plotted against Y « 2 :Yt . ii for a fixed time T 
for a number of coordinate values X, and the position of the 
values >£ a constant read from these curves. In figure 12, the 

cylinder is shown from X - jt/2 to X » it. The separation time is 

given by 2T 2 00 s X « -2.3 1 !, that is, the start of the separation 
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by T = I.O 80 In figure 12 2T^ = 5 j lienee T - 1 . 58 , "was chosen. 

For this chosen time, the separation point has already progressed up 
to beyond 60 ° at the cylinder; nevertheless the boundary layer still 


is fairly thin, the relative sizes correspond to the values S = 10 cm, 

2 

k = 0 . 0 l£S_ (water), Y = 0.1 PS —, that is, to a very small acceleration, 
sec a -„2 


Accordingly, t = 15.8 sec. 




Figure 12 
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The picture obtained by the 

for V * 10 - 2 S- after 1 .58 sec. 
sec* 1 

thickening of the boundary layer 
of 1 :^L 0 . 


previous reduction formulas 
is represented in figure 12 . The 
vould be diminished in the ratio 
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